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Abstract 


Precipitation is the main component of the hydrological cycle, which is very important in many 
studies, including those focusing on management of irrigation and drainage systems, crop 
performance, environmental issues, floods, droughts, etc. Therefore, this study eveluated seven 
GCMs models of Coupled Model Intercomparison Project phase 6 (CMIP6) regarding the 
simulation and projection of precipitation changes in the northwest of Iran under three scenarios of 
SSP in three periods (2021-2050, 2051-2080, 2081-2100) and compared them to the historical 
period (2014-1985). The trend of precipitation changes was calculated using Mann-Kendall test 
and Sen’s slope estimator. The Linear Scaling Bias Correction (LSBC) was used to downscaling 
the GCMs data and the validity of different models was evaluated using RMSE, MAE and R2 
indices. The results of performance evaluation of CMIP6 models showed that among the studied 
models, MPI-ESM1-2-LR model with an average R2 of 0.86 and RMSE equal to 19.7 at the 
regional level were more accurate than other models in the simulation of precipitation. The results 
of projection also showed that precipitation according to SSP1-2.6 scenario will increase by 2.6% 
in the next three periods and on average in the end of 21st century and at most stations it is 
significant at the level of 0.01%. But according to SSP3-7.0 and SSP5-8.5 scenarios, the amount of 
precipitation will decrease by 14.5% and 3.6% in the end of 21st century, respectively, with the 
highest decrease being related to the rainy areas in the southwest of the region. In general, 
according to different scenarios, in most of the studied areas, the precipitation trend is decreasing 
until the end of 21st century, and it is necessary for officials and planners to adopt the required 
strategies to adapt to the resulting climate change. 
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حکیده 


بارش اصلی‌ترین جزء چرخه هیدرولوژی است که در بسیاری از مطالعات ازجمله 
مديريت سیستم‌های آبیاری و زهکشی. عملکرد محصولات کشاورزی. مطالعات 
زیست محیطی. سیلاب. خشکسالی و غيره از اهميت زيادى برخوردار است. لذا در اين 
پژوهش به ارزیابی Y‏ مدل GCMs‏ از CMIP6‏ در شبیه‌سازی بارش و پیش‌نگری 
تغییرات بارش در شمال غرب ایران تحت سه سناریوی واداشت اجتماعی-اقتصادی 
(SSP)‏ در سه دوره (۲۰۲۱-۲۰۵۰ ۲۰۵۱-۲۰۸۰ ۲۰۸۱-۲۱۰۰) نسبت به دوره پایه 
(۱۹۸۵-۲۰۱۶) پرداخته شد. روند تغییرات بارش با استفاده از آزمون من کندال و 
تخمینگر شیب سن محاسبه گردید. جهت مقیاس کاهی داده‌های 5 نیز از رویکرد 
مقیاس بندی خطى (LSBC)‏ استفاده شد و صحت سنجی مدل‌های مختلف با استفاده 
از سنجه‌های متداول آماری» MAE RMSE‏ و 82 ارزیابی گردید. نتایج حاصل از 
ارزیابی عملکرد مدل‌های CMIP6‏ نشان داد که از ميان مدل‌های موردبررسی» مدل 
MPI-ESMI2-LR‏ با متوسط R2‏ برابر با ۰/۸۲ و RMSE‏ برابر با ۱۹/۷ در سطح 
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نتایج حاصل از پیش‌نگری تغییرات بارش نيز نشان داد بارش بر اساس سناریوی 
SSP1-2.6‏ در هر سه دوره odul‏ و به‌طور متوسط تا پایان قرن ۲۱ به میزان 7/1 درصد 
افزایش بيدا خواهد کرد و در بیشتر ایستگاه‌ها در سطح ۱ درصد معنی‌دار شده 
است؛ اما بر اساس سناریوهای SSP3-7.0‏ و SSP5-8.5‏ میزان بارش تا پایان قرن ۲۱ به 
ترتیب ۱8/۵ و ۳/۰ درصد کاهش lay‏ خواهد کرد که بیشترین میزان کاهش أن نيز 
مربوط به مناطق پر بارش در جنوب غرب منطقه است. میزان کاهش در مناطق غربی و 
جنوت غربى در سطح v/v)‏ درصد معنی‌دار شده As]‏ درمجموع با توجه به 
است و لازم است مسئولان و برنامه ریزان راهکارهای لازم جهت سازگاری با تغییرات 
اقلیمی حاصل را اتخاذ نمایند. 


کلیدواژه‌ها: بارش, تغيير اقليم» شمال غرب. مقیاس بندی خطی. CMIPOSSPS‏ 


١‏ - مقدمه 


كرمايش جهانى در حال حاضر يذيرفته شده است (ژانگ و همکاران » ۲۰۱۷). یکی از نمودهای كرمايش 
جهانى و تغيير اقليم افزايش فراوانى و شدت دورههاى همراه با دماهاى بالا و يايين است و برخى مطالعات 
نشان دادهاند كه تأثيرات فرينهاى اقليمى آشکارتر و مستقيمتر از ميانكينهاى آب‌وهوایی است (هيأت بين- 
الدول تغيير اقلیم » ۲۰۱۳). دمای بالاتر نگهداشت رطوبت جو را افزايش داده و اين امر به نوبه خود می تواند 
stabil ay ona‏ بارشو فاق کین قوف OL)‏ و OHA Ob aa‏ فك انيف DIR‏ کات و ting‏ 
فضايى رويدادهاى فرين آب‌وهوایی افزايش LL‏ (هورتون و همکاران ۲۰۱۵). مسئله تغيير اقليم امروزه 
به‌عنوان یکی از شايعترين مباحث علمى و حتى سياسى-اجتماعى مطرح است. اصولاً تغيير و حركت جزء 
وجودى سامانه‌های طبيعى است. انسان در تلاش است كه با ابزارها و دانستنىهاى خود قانونمندی اين 
تغييرات را کشف كند و بدينوسيله اقدام به پیش‌بینی رفتارهای بعدى سامانه آبوهوا نمايد (عزیری,۱۳۸۳). 
فرايند تغییر اقليم بهويزه تغييرات دما و بارش مهمترين بحث مطرح در قلمرو علوم محيطى می‌باشد. اين يديده 
به دلیل ابعاد علمى و كاربردى (اثرات محيطى و اقتصادى- اجتماعى) آن از اهميت فزاینده‌ای برخوردار است» 
جراكه سيستمهاى انسانى وابسته به عناصر اقليمى مانند کشاورزی. صنايع و امثال آن بر مبناى ثبات و پایداری 
اقلیم طراحى شده و عمل می‌نمایند (ذهبيون و همكاران. (YA‏ در اين بين نيز بارش مهمترين متغير 
Zhang et al‏ 1 

2 IPCC 


3 Yang et al 
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سال سيزدهم پیش نكرى بارش با استفاده از مدلهاى 081186 تا يايان قرن 7١‏ در dius‏ غرب ايران Wo‏ 


آب‌شناختی است که ييوند ميان جو و فرآیندهای سطحی را برقرار می‌سازد و ناهنجاری‌های مثبت و منفى آن 
سبب بروز سيل و خشكسالى می‌شود (مسعوديان و همكاران. ۱۳۹۳). برای ارزيابى تغييرات اقليمى 
مدل‌های جهانى اقليم (GCMs)‏ طراحى شدهاند. اين مدل‌ها خصوصيات فیزیکی. چرخش‌ها و حركات 
اتمسفرى را تحليل مىنمايند و به دنبال آن متغيرهاى هواشناسی را در شبکه‌های خاص شبيهسازى می HES‏ 
همجنين اين مدلها در حال حاضر قوىترين ابزار به‌منظور توليد سناريوهاى اقليمى می‌باشند (شائمى و 
حبیبی نوخندان. (AYM‏ مدل‌های GCM‏ متعددى توسعه يافته است و به دليل تفاوت ساختار اين مدلها و 
تفاوت شرايط اولیه‌ای که تحت آن اجرا می‌شوند. نتايج آنها حتى برای یک سناريوى انتشار OLS‏ متفاوت 
است. لذا برای استفاده در مطالعات منطقه‌ای» لازم است خروجی اين مدلها با استفاده از یک روش مناسب 
ريز مقياس شود (كى و همكاران » 9 با توجه به اهميت بارش لازم است اين يارامتر در مناطق مختلف 
برآورد شود. پیش‌بینی بارش برای مناطقى كه اقتصاد آن‌ها مبتنى بر كشاورزى است» بسيار مهم است. از 
طرفى با بيش بينى بهموقع بارندكى امكان مقابله با سیلاب» خشک‌سالی و كاهش خسارت ناشی از آن‌ها نيز 
مهيا می‌شود. در اين راستا فريرا و همكاران' (۲۰۱۸) به پیش‌نگری تغييرات بار شهاى تابستانى در جنوب 
شرقى ابالات‌متحده پرداختند. نتايج حاكى از افزايش قابل‌ملاحظه بارش به دليل افزايش دما و شار بخارآب 
وف SS‏ فیکار (P989‏ وی الراك کیرات Clays DU de Los‏ پو تا در عند با Jald‏ 
سناريوهاى RCP‏ پرداختند. نتايج OUS‏ داد دما و بارش در دوره آینده افزايش بيدا خواهد كرد. جیانگ و 
همکاران* (۲۰۲۰) به پیش‌نگری تغييرات بارش در آسياى مركزى بر اساس ۱۵ مدل GCM‏ از “CMIP6‏ تحت 
٤‏ سناریوی ۹5۳" پرداختند. نتايج نشان داد بر اساس تمام سناريوها ميانكين بارندگی سالانه افزايش خواهد 
پافت. ژانگ و همكاران" (۲۰۲۱) قابليتهاى ۶۲ مدل از CMIP6‏ را جهت بازتوليد دما و بارش در سين 
کیانگ جين را مورد ارزيابى قرار دادند. نتايج OUS‏ داد اين مدلها توانايى خوبى در شبيهسازى دما و بارش 
دارند و استفاده از cus‏ جند مدل نيز عملكرد بهترى نسبت به مدلهاى منفرد دارند. مجدى و همکاران" 
(۲۰۲۲) به پیش‌نگری تغييرات دما و بارش در خاورميانه و شمال آفریقا با استفاده از متوسط ۲۳ مدل GCM‏ 
و دو سناريوى SSP‏ يرداختند. نتایج SE‏ از افزايش دما و كاهش بارش در بيشتر مناطق موردمطالعه است. 


صلاحی و همکاران (۱۳۹۵) به پیش‌بینی تغييرات دما و بارش با استفاده از داده‌های خروجى مدل 
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3 تحت سه سناريوى انتشار در حوضه آبريز درياجه پرداختند. نتایج حاكى از کاهش ۸/۲ درصدی 
بارش در دهه ۲۰۵۰ در حوضه موردمطالعه است. خورشيددوست و همكاران (۱۳۹۷) به پیش‌بینی بارش هاى 
سنگین در سواحل جنوبی دریای خزر در دوره ۱۳۸۹ تا ۱۶۰۸ پرداختند. quU‏ نشان داد در دوره آتی 
درمجموع بر تعداد روزهای همراه با بارش سنگین در اکثر ایستگاه‌های موردمطالعه افزوده خواهد شد. خزانی 
و همکاران (۱۳۹۸) با استفاده از مدل HadGEM2‏ به پیش‌نگری تغییرات دما و بارش در حوضه سیرا 
پرداختند. نتایج نشان داد كه در دوره آینده بارش کاهش و دما افزایش می‌یابد. حیدری و همکاران(۱۳۹۸) 
چشم‌انداز تغییرات بارش را در منطقه شمال غرب با استفاده از مدل‌های گردش عمومی جو و سناریوهای 
RCP‏ مورد ارزیابی قرار دادند. نتایج نشان داده که بیشتر مدل‌ها برای دوره‌های آینده افزایش بارش را پیش - 
بینی می‌کنند. کثیری و همکاران (VE)‏ به پیش‌نگری تغییرات بارش در سواحل جنوبی دریای خزر با 
استفاده از سناریوهای مختلف واداشت تابشی پرداختند. نتایج نشان داد بارش در سواحل جنوبی دریای خزر 
تا پایان قرن ۲۱ نسبت به دورة پایه به‌طور متوسط ۵/۷ درصد افزايش خواهد یافت. 

جمع‌بندی پیشینه تحقیق نشان داد که بیشتر مطالعات انجام شده در زمینه تغییر اقلیم در ایران و در منطقه 
شمال غرب با استفاده از برونداد مدل‌های گزارش سوم (CMIP3)‏ و يا گزارش پنجم (CMIPS)‏ و داده‌های 
اين مدل‌ها صورت گرفته است و بررسی بارش به‌صورت جامع و با استفاده از مدل‌های گزارش ششم 
6 صورت نگرفته است. لذا در اين پژوهش درستی سنجی مدل‌های متداول و پرکاربرد گزارش ششم 
در دسترس با استفاده از متداول‌ترین سنجه‌های آماری بررسی شد و تغییرات بارش در شمال غرب ايران در 
دوره آینده (۲۰۲۱-۲۱۰۰) با استفاده از خروجی مدل منتخب با کمترین میزان خطا تحت سناریوهای مختلف 
SSP‏ بيش نكرى گردید. 

با توجه به آن که منابع آب در معرض خطرات ناشی از تغييرات اقلیم قرار دارند. بررسی تغییرات بارش در 
سال‌های آینده می‌تواند راهگشای معضلاتی همچون خشکسالی. سیلاب‌های ناگهانی و تخریب‌های 
زیست محیطی باشد. 


Y‏ - داده‌ها و روش‌شناسی 
۱-۲- منطقه موردمطالعه 


محدوده موردمطالعه اين پژوهش. شمال غرب ol ul‏ است که از عرض جغرافیایی ۳۶ درجه و ٤٤‏ دقیقه 
الى YA‏ درجه و ۶۷ دقيقه عرض شمالی و طول جغرافیایی ٤٤‏ درجه و ۵ دقیقه الى £A‏ درجه و ۵۲ دقیقه 


داده‌های بارش Yo‏ ايستكاه همديدى به‌صورت روزانه» با توزيع مناسب و دارای دوره آماری بلندمدت 


سال سیزدهم پیش نگری بارش با استفاده از مدل‌های 0811276 تا SLL‏ قرن ۲۱ در شمال غرب ایران ۱۷ 


(۱۶ ۱۹۸۵-۵۰ حداقل سی‌ساله در سطح منطقه موردمطالعه استفاده LŚ‏ كه موقعيت جغرافيايى Lal‏ در 
شکل (۱) آورده شده است. 
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شكل -١‏ موقعيت جغرافیایی و توي وكرافى منطقه موردمطالعه و ایستگاه‌های هواشناسی موردبررسی 


۳-۲- مدل جهانى (GCM)‏ و سناريوهاى موردبررسى 

در اين پژوهش از دادههاى بارش ایستگاه‌های همديد شمال غرب كشور به‌صورت روزانه در دوره آماری 
تلد و ساسا با کرو قاری ستل‌های جما 548-9313 اقات كلد که ان ماران هر اکتا 
کشور دریافت گردید. جهت بررسی تغییرات بارش در دوره‌های آینده موردبررسی (۲۰۲۱-۲۱۰۰) ابتدا ۷ 
مدل GCM‏ از CMIP6‏ که دارای داده‌های بارش با وضوحبالا بودند اتتخاب شدند. مشخصات مدل‌های 
موردبروسی در جدول (۱) آوزده شده است, داده‌های بارش مدل‌های مذ‌کور از مرکر ۳50 دانلود گردید, 
سپس polie‏ بارش مشاهداتی و تاریخی برای هر یک از ایستگاه‌های موردمطالعه با تهیه پرنامه در محيط 
0 با استفاده از نزدیک‌ترین كريد داده‌های GCMs‏ در دوره پایه (۱۹۸۵-۲۰۱۶) استخراج 
گردید. سپس از روش تصحيح اریبی يا روش مقیاس کاهی خطى جهت ریزگردانی داده‌ها استفاده گردید. در 
ا ی sales‏ و RUN E‏ ا ای رر اکن ریات غا 


1 https://esgf-node.lInl.gov/search/cmip6/ 


\VA‏ جغرافيا و مخاطرات محيطى شمارة اول 


(RMSE)‏ ميانكين مطلق خطا (MAE)‏ و ضريب تعیین Rsqrt)‏ ارزيابى كرديد. يس از ارزيابى مدلهاى 
مختلف GCM‏ با استفاده از بهترين مدل» تغييرات بارش در دوره آينده (۲۰۲۱-۲۰۵۰) تحت سه سناريوى 
واداشت اجتماعى- اقتصادی (SSP)‏ پیش‌نگری كرديد و تغييرات آن در دوره آينده نسبت به دوره پایه 
محاسبه گردید. دليل انتخاب مدل‌های مذكور از مجموعه مدلهاى CMIP6‏ اعتبار دقت و در دسترس بودن 
سناریوهای موردنیاز جهت مطالعه است. ينج سناریوی SSP‏ از 5581 تا 85۳5 وجود دارد که به ترتیب 
نشان‌دهنده گذر از جاده سبز» وسط جاده جاده سنگی. جاده تقسیم شده و گرفتن بزرگراه است. کمیت‌های 
SSP‏ اساس تلاش مشترک بين IAV‏ و جامعه IAM‏ ساخته شده است. سناریوهای موردمطالعه در اين 
پژوهش يعنى SSP126‏ ۹5۳265 و 55۳585 نشان‌دهنده انتهای يايين محدوده مسیرهای اجباری آینده 
انتهای متوسط تا بالا محدوده مسیرهای اجباری آينده و انتهای بالای محدوده آینده هستند (زهرائی و 
حسینی, ۱۳۹۹). در ادامه با استفاده از آزمون آماری من کندال و تخمینگر شیب سن, روند و شيب تغییرات 
بارش در منطقه شمال غرب کشور مشخص گردید. 


Pn‏ مدل موسسه/ کشور قدرت تفکیک 

MRI-ESM2-0 ۱‏ موسسه تحقيقات جوی» ژاپن APT‏ 

۰/4/4 موسسه ماكس پلانک. آلمان‎ MPI-ESM1-2-LR Y 
Yox) آزمایشگاه زئوفيزيك و دینامیک سیالات ایالات‌متحده‎ GFDL-ESM4 ۳ 
۱۰۶ ۱ oe el MIROC6 t 
۷۲/۸۸ و‎ M CL d MIROC-ES2L 0 
۱۱ جين‎ om موسسه فیزیک‎ FGOALS-g3 1 

CanESM5 ۷‏ مرکز تحلیل و مدل‌سازی اقلیمی. کانادا ۱/۱۱ 


منبع: ژانگ و همکاران (۲۰۲۱) 
-tY‏ مقیاس AS‏ 


جهت مقیاس کاهی داده GCMs‏ از روش مقیاس GAS‏ خطى استفاده گردید. مقیاس کاهی بر اساس 
ميانكين تفاوت بين سری‌های زمانى مشاهده شده ماهانه و سری‌های زمانی دوره تاریخی GCM/RCM‏ در 
همان دوره از سرى زمانى مشاهده شده است. سپس اين تفاوت‌ها براى داده‌های آب‌وهوایی شبیه‌سازی شده 


برای دریافت متغیرهای آب‌وهوایی با اریبی اصلاح شده اعمال می شوند (شرستا و همکاران۲۰۱۵). رويكرد 


1 root-mean-square error 
2 mean absolute error 


سال سيزدهم پیش نكرى بارش با استفاده از مدلهاى CMIPO‏ تا يايان قرن 7١‏ در شمال غرب ايران WA‏ 


مقياس بندی خطی (LSBC)‏ با مقادير اصلاح شده ماهانه بر اساس تفاوت بين مقادير مشاهداتى و 
شبیه‌سازی شده عمل مىكند. داده‌های بارش با ضریبی بر اساس نسبت ميانكين بلندمدت ماهانه داده‌های 
مشاهداتى و داده‌های کنترل اصلاح مى شود. در op!‏ حالت شرط بر اين است که ضرایب تصحیح اعمال شده 


وافزودهها حتى برای شرايط آتى نيز بدون تغيير باقی می‌مانند (توجباين و AYN Ole‏ 


Um Pops (t) \ 
Pee USP ee cy, | مستت‎ | Q) 


Hm Pcontr (t) 


PRE) = Po [East 00 


Hm Pcontr (t) 


در روابط P YL‏ بارش است؛ contr‏ تعداد سرىهاى شبیه‌سازی شده GCM‏ در طول »55« کنترل Obs‏ 
دوره زمانى مشاهداتى fre‏ سرس زمانى پیش‌نگری شده آينده است كه بايد اريبى ol‏ حذف شود؛ BC‏ سرى 
زمانى پیش‌نگری شده آينده است كه اريبى آن حذف شده است. t‏ گام زمانى و Hm‏ ميانكين بلندمدت ماهانه 
است (منذر و (Y Yr col s‏ با اعمال روش مقياس كاهى بر روى دادههاى شبیه‌سازی شده عملکرد مدل 


در شبيهسازى بلندمدت متغير بارش به مقدار زيادى افزايش مىيابد. 
؟-0- آزمون نايارامتريك آمارى من کندال 


آزمون نايارامترى من- كندال که توسط من در سال ١950‏ ارائه و سپس توسط كندال در سال ۱۹۷۵ 
تكميل گردید. بر Gal‏ مرتبه داده‌ها در یک سرى زمانی استوار است. اين آزمون برای بررسى تصادفى بودن 
دادهها (عدم وجود روند) در مقابل وجود روند در سرىهاى زمانى هيدرولوزيكى و هواشناسى مورد استفاده 
قرار م ى كيرد (ژانگ و همکاران. ۲۰۰۰). Cope‏ اين آزمون نسبت به ple‏ آزمون‌های تعيين روند استفاده از 
مرتبه داده‌ها در سرى زمانى بدون در نظر داشتن مقدار متغيرها می‌باشد كه به دليل وجود جنين خاصيتى» 
می‌توان از اين آزمون برای داده‌های داراى جولكى نیز استفاده كرد و نيازى نيست كه دادهها در قالب توزيع 
خاصى در آیند. فرض موردبررسی در اين آزمون به‌صورت زیر می‌باشد: 

Ho‏ داده‌ها به‌صورت تصادفی توزیع شده‌اند (روند ندارند). 

Hy‏ داده‌ها دارای روند هستند. 


در اين i gal‏ اگر X1,X2,.--5Xn‏ مشاهدات مورد نظر باشند» آنگاه داریم: 


1 Linear Scaling Bias Correction 


۱۸۰ جغرافيا و مخاطرات محيطى شمارة اول 


71-1 n 


uS Y'sm(s, - x,) (v) 


kal j=k+1 
تابع علامت می‌باشد كه به‌صورت زیر تعريف می‌گردد:‎ sgn كه رابطه فوق 1" مقدار داده زام و‎ 


+1 if (x, -x,)»0 
sgn(x)- 40 if (x, -x,)-0 (£) 
-1 dr (x, - x, )«0 
در واقع در این آزمون هر داده با تمامی داده‌های پس از خود مقایسه می‌شود؛ در این مرحله مىتوان‎ 
بجای استفاده از مقادير اصلی دادههاء از مرتبه داده‌ها درمجموعه مورد نظر (سری زمانی) استفاده کرده و‎ 
(Rank base) مرتبه‌ها را به همین روش مقایسه نمود. به اين دلیل اين آزمون یک آزمون بر پایه مرتبه داده‌ها‎ 
می‌باشد. با فرض این که داده‌ها مستقل بوده و توزيع يكنواخت دارند. ميانكين و واریانس 5 از روابط زیر به‎ 


دست مىآيد: 


E(S) =0 (0) 


m 


n(n -1(2n 4 5)- X r(t, -1)(2t, +5) 5 


Var(S)= T 


در رابطه بالا: 

n‏ تعداد دادههاء m‏ تعداد گره‌ها و t‏ تعداد osla‏ در هر گره می‌باشد. 

منظور از گره اين است که اگر از یک مقدار داده, بیشتر از یکی وجود داشته باشد. این polite‏ مساوی. 
تشکیل یک گره را می‌دهند و تعداد اين polio‏ مساوی در ply elm oS‏ ) می‌باشد. 


آماره اين آزمون (Z)‏ داراى توزيع نرمال بوده و از رابطه زیر به دست می‌آید: 


i S>0 
JVar(s) f 
Z=40 if S=0 (V) 
S+1 s<0 


فرض صفر يذيرفته مىشود و در غير اين صورت. فرض صفر رد خواهد شد. در حالت رد فرض صفر 


سال سيزدهم پیش نكرى بارش با استفاده از مدلهاى 081156 تا يايان قرن 7١‏ در dius‏ غرب ايران ۱۸۱ 


(وجود روند)» درصورتی که 0 باشد. سری زمانی دارای روند مثبت (صعودی) و درصورتی که 0 باشد. 
|Z| > 1.64‏ در سطح معنی‌داری E‏ درصد. اگر 


zo Z| 2.57 <€‏ . 
در سطح معنی‌داری ۵ درصد و اگر pe‏ ادر سطح معنىدارى ۱ درصد فرض صفر (عدم 


وجود روند) رد مى شود A(\AV0 JS)‏ 


سری زمانى داراى روند منفى (نزولى) خواهد بود. اگر 
|Z| > 1.96‏ 


esas 
برای تشخيص روند در یک سرى زمانى داده‌ها ارائه شد. اين‎ (OAM) سن‎ dau تخمين شيب خط روند‎ 
روش یک روش نايارامترى است که با استفاده از تحليل تفاوت بين مشاهدهاى یک سرى زمانی» یک شيب‎ 
را برای سری زمانی داده‌ها از رابطه زیر محاسبه می‌کند. در روش سن برای هر دو جفت داده متوالی‎ the 
شیب خط‎ laol یک شیب محاسبه می‌شود. شیب‌های محاسبه شده تشکیل یک سری زمانی می‌دهند که میانه‎ 
می‌دهد. در آزمون من- کندال شیب خط روند سری داده‌ها با روش تخمین گر سن از رابطه‎ OLY روند را‎ 


ناپارامتری زیر محاسبه می‌گردد. 


"UN ar c PP 
B = Median | ۳ > i) 
J-1 (A) 

كه در B ol‏ برآوردگر شيب خط روند و زک به ترتيب polis‏ مشاهداتى 1ام j‏ ام می‌باشند و ريك 


واحد زمانى بعد از i‏ می‌باشد مقادير مثبت و منفى 8 نشان دهنده روند افزايشى و كاهشى در سرى داده- 
هاست (سن,» 8 


۷-۲- ارزيابى و صحت سنجى عملكرد مدل‌ها 

يس از استخراج خروجى مدلها براى صحت سنجی آن‌ها از سنجه‌های متداول آمارى استفاده شد كه در 
ادامه شرح مختصرى از اين شاخصها آورده شده است. در رابطه (A)‏ كه تحت عنوان ضريب تعیین (Rsqtt)‏ 
شناخته مى شود معيارى بدون بعد است و بهترين مقدار آن برابر يك می‌باشد (صداقت کردار و همکاران 
(WV ۷‏ مجذور ميانكين مربعات خطا (RMSE)‏ به‌عنوان قياسى برای OUS‏ دادن اختلاف بين مقادير 
شبیه‌سازی شده از مقادير اندازه‌گیری به كار مىرود اين معيار كه به‌صورت رابطه (۱۰) تعريف می‌شود 


به‌عنوان مرسومترين شاخص خطا به كار مىرود (لين و همکاران» 2007. ميانكين مطلق خطا (MAE)‏ برای 


1 Sen's Slope Estimator 


\AY‏ جغرافيا و مخاطرات محيطى شمارة اول 


مقایسه عبارت به عبارت خطای نسبی مقادیر شبیه‌سازی شده با توجه به مقادير اندازه‌گیری شده به كار مى- 


رود كه بە‌صورت رابطه (۱۱) ارائه م $ ,23 (هو و همكاران» Que‏ 


00 
R? 55 i=l 
N 
20, 2K, 
5 (4) 
FN 
> مع‎ - X.Y 
RMSE = (| -= 
N (V+) 
N 
۳2 9X 
MAE = id 
N QD 


در روابط فوق Xo‏ داده‌های مشاهداتى» X,‏ داده‌های شبيهسازى شده» Xo‏ ميانكين داده‌های مشاهداتىء. N‏ 


تعداد داده‌ها می‌باشد. 


۳- نتایج و بحث 


ay gut‏ دست آمده از اعمال آزمون من-کندال و تخعمینگر شیب سن بر روی داده‌های بارش سالائه در 
دوره آماری GL‏ (۱۹۸۵-۲۰۱۶) در ایستگاه‌های موردمطالعه نشان می‌دهد که بارش در ایستگاه‌های اردبیل» 
يارس آباد. خوی» جلفاء زنجان» سراب. سردشت و ماكو دارای روند افزایشی و در ساير ایستگاه‌ها دارای روند 
کاهشی است. بر اساس آزمون من کندال روند افزایشی بارش در هيج كدام از ایستگاه‌های موردمطالعه معنی- 
دار نشده است و تغییرات مذکور از نوع نوسانات کوتاه‌مدت آب‌وهوایی در منطقه می‌باشد؛ اما روند کاهشی 
در ایستگاه‌های سنندج و مراغه در سطح ۰/۰۱ درصد و در ایستگاه‌های سهند و مریوان در سطح ۰/۰۵ درصد 
معنی دار شده است و حاکی از تغييرات کاهشی بارش در ایستگاه‌های مذكور است. بیشترین شیب تغییرات 
کاهشی نيز مربوط به ایستگاه مریوان با شیب سن برابر ۱۹/۲- است. بیشترین شیب تغییرات افزایشی نيز 
مربوط به ایستگاه سردشت با 1/۵ است. درمجموع تغییرات کاهشی بارش در منطقه موردمطالعه بیشتر از 


تغییرات افزایشی است. (جدول AY‏ 
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جدول ۲- نتايج حاصل از آزمون‌های آمارى بر روى سرىهاى بارش سالانه ایستگاه‌های موردمطالعه 


رديف ايستكاه ضريب من كندال شيب سن رديف ايستكاه ضريب من كندال شيب سن 
١‏ اردبيل M YS WAR!‏ زرینه ۱/۹۲ ۳۱/- 
۲ ارومیه ۰/۲۹۰ -£6/* IE‏ زنجان ۰/۵۶ ۰/۸۵ 
۳ اهر M S VM-‏ سراب YAV Vit‏ 
3 بيجار +/4V-‏ -£0/\ 1۷ سردشت UNY‏ 1/۵ 
۵ يارس آباد JAY‏ ۱/۳۳ ۸ سقز ۱/۸ 4/47- 
1 پیرانشهر ۰/۲۸۰ ۰/۷۲ ۱۹ سنندج -—VYA HRY PY‏ 
۷ تبریز AD‏ ۱/۲۸۰ ۳۰ سهند -۲/۰۳۲ # Y/* Y4-‏ 
۸ تکاب ۱/۳۷ ۲/۸۹ ۲۱ قروه -۰/۷۵ -MYY‏ 
4 جلفا 1/0۸ ۳/۳۰ ۳۲ ماكو A‏ ۳/۹ 
E‏ خرمدره اكه UYM‏ ۲۳ مراغه -0/V4 RY /AQ-‏ 
Mw‏ خلخال -۱/۸۵ -۱/۵۲ yé‏ مريوان *YY-‏ ۷- 
ZA iY Ce Ww‏ 
Yo‏ مهاباد -۰/۹۳ -Y‏ 
۱۳ ميانه -14/+ ZA‏ 
#معنی‌داری در سطح ۱ درصد #معنی‌داری در سطح 6 درصد 


توزیع فضایی تغییرات بارش نيز بر اساس آزمون من کندال و تخمینگر شیب سن نشان می‌دهد که به‌جز 
oly‏ شمالی منطقه مورد و منطقه ارسباران و همچنین مناطقی به‌صورت پراکنده در شرق و غرب منطقه در 
yhe‏ مناطق بارش دارای روند کاهشی است که بیشترین تغییرات أن نيز مربوط به نواحی جنوبی و جنوب 
غربی منطقه موردمطالعه و همچنین جنوب شرق دریاچه ارومیه است. 

در اين پژوهش, ۷ مدل GCMs‏ مورد ارزیابی قرار گرفت و پس از اصلاح اریبی به روش 1.5130 به 
cli!‏ عملگرد مد ها يد بروسی etal‏ بارش Bilt, Sala,‏ مدلها در در ا (Y8483Y-15)‏ 
پرداخته شد. نمونه‌ای از اصلاح اريبى بر روى دادههاى بارش ايستكاه تبريز در شکل (۳) آورده شده است. 
همانطور كه مشاهده مىشود دادههاى بارش داراى عدم قطعيت زيادى است لذا انجام عمل اصلاح اريبى 
توانسته است به ميزان زيادى عملكرد مدل را در شبيهسازى بلندمدت بارش در منطقه موردمطالعه افزايش 
دهد. با توجه به تعداد زياد ایستگاه‌ها و نمودارهاى حاصل» فقط خروجى ايستكاه تبريز به‌صورت نمودار در 
اين بخش ارائه كرديد. برای مطالعه بيشتر مىتوان به مطالعات شرستا' (۲۰۱۵) و تتجبين' (Y)‏ در اين 


زمينه مراجعه كرد. 


1 Shrestha 
2 Teutschbein 
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— Corrected 


شکا polis -Y‏ مشاهداتی. تاریخی و اصلاح شده بارش در ايستكاه تبريز در دوره ak‏ (۱۹۸۵-۲۰۱۶) 


يس از مقياس كاهى به نتايج حاصل از ارزيابى مدلهاى مختلف با استفاده از شاخص مختلف در 
ایستگاه‌های موردمطالعه پرداخته شد. نتايج نشان داد دقت اكثر مدل‌های موردبررسى در مناطق كم بارش و 
پربارش منطقه موردمطالعه كمتر از مناطق دیگر است و می‌توان بیان كرد تغييرات زياد بارش توي وكرافى و 
سيستمهاى جوى مختلف از مهم ترين دلايل افزايش خطا در برخى ایستگاه‌ها و مناطق موردمطالعه است. با 
توانند الكوى بارش را به‌طور دقيق شبیه‌سازی كنند. درمجموع و بر اساس R? polis‏ مدل MPI-ESM1-2-LR‏ 
با داده‌های بارش مشاهداتی در منطقه موردمطالعه برابر با ۰/۸۹ و ميزان RMSE‏ أن نيز برابر با ۱۹/۷ است. 
نتایج رتبه‌بندی عملکرد مدل‌های موردبررسی در شبیه‌سازی بارش شمال غرب کشور در جدول Y‏ آورده شده 
است. يس از مدل MPI-ESMI-2-LR‏ دو مدل MIROC-ES2L‏ و MIROC6‏ هر دو با ضریب تعیین برابر با 
۵ بهترین مدل‌های CMIP6‏ برای پیش‌نگری بارش در منطقه موردمطالعه هستند. همچنین مدل‌های 
MIROC-ES2L‏ و FGOALS-g3‏ در حد متوسط هستند. مدل CanESMS‏ نيز با R?‏ برابر با ۰/۶۹ ضعیف- 
ترين عملكرد را در بين مدلهاى موردبررسى جهت شبيهسازى بارش در منطقه موردمطالعه داشته AL‏ بر 
اساس شاخحص RMSE‏ و MAE‏ نيز مدل MPEESMI-2-LR‏ بهترين عملكرد را در بين مدلهاى 
موردبررسی داشته aul‏ ولى در برخى مناطق دارای کم برآوردی و در برحی مناطق pe‏ دارای بیش 
برآوردی بوده است. بر اساس دو شاخص مذکور نیز مدل FGOALS-g3‏ با RMSE‏ برابر با ۶۰/۳ و MAE‏ 


سال سيزدهم پیش نكرى بارش با استفاده از مدلهاى CMIPO‏ تا يايان قرن 7١‏ در شمال غرب ايران ۱۸۵ 


جدول ۳- رتبه‌بندی مدل‌های 021176 جهت شبيهسازى بارش بر اساس ميانكين سنجه‌های خطاسنجی 


ردیف مدل ضريب تعيين | مجذور میانگین مربعات خطا | میانگین مطلق خطا 
MPI-ESM1-2-LR ۱‏ كلل ۱۹/۷۱ ۱9/۰ 
MAE YY/VO */A0 MIROC-ES2L Y‏ 
YY/AA ۳۸/۰۵ ۰/۸۵ MIROC6 ۳‏ 
VINE AY MRI-ESM2-0 t‏ ۱۹/۸۹۶ 
YYAA ۰/۲ (MO FGOALS-g3 ۵‏ 
YAAA YVOA "oY GFDL-ESM4 1‏ 
CanESMS M‏ ۹ 10/4۱ ۲۰/۳۹ 


يس از انتخاب بهترين مدل با بهترين عملكرد در شبیه‌سازی بارش منطقه موردمطالعه يعنى مدل MPI-‏ 
ESMI-2-LR‏ به پیش‌نگری بارش با استفاده از مدل مذكور در دوره‌های آینده (۰۲۰۲۱-۲۰۵۰ ۲۰۵۱-۲۰۸۰ 
و ۲۰۸۱-۲۱۰۰) براساس سه سناريوى موردبررسى يرداخته شد. نتايج نشان داد كه بر اساس سناریوی 
5581-6 بارش در منطقه موردمطالعه در دوره (۲۰۲۱-۲۰۵۰) بين ۲۰۰ تا ۸۵۰ میلی‌متر خواهد بود كه 
نسبت به دوره يايه به‌طور متوسط ۳/۲ درصد YN)‏ ميلىمتر) در منطقه افزايش OUS‏ مىدهد. در دو دوره 
دیگر نيز (۲۰۵۱-۲۰۸۰ و ۲۰۸۱-۲۱۰۰) ميزان بارش بهطور متوسط به ترتيب ٠/05‏ و ٤/٦‏ درصد نسبت به 
دوره يايه افزايش خواهد يافت كه نسبت به دوره قبلى كاهش زيادى را نشان مىدهد. بر اساس توزيع فضايى 
بارش در دوره‌های آينده در هر سه دوره موردبررسى بيشترين تغييرات افزايشى مربوط به بخش‌هایی شمالى 
منطقه موردمطالعه يعنى منطقه جنگلی‌های ارسباران است كه می‌تواند به دليل يوشش جنگلی منطقه» نزديكى 
به منابع آبی و همجنين توي وكرافى منطقه باشد. در ساير مناطق موردمطالعه ميزان بارش كاهش پیدا خواهد 
كرد كه بيشترين ميزان آن مربوط به نواحى جنوب غربى منطقه موردمطالعه مىباشد كه جزو مناطق يربارش 
در منطقه موردمطالعه نيز هستند. درمجموع ميزان بارش بر اساس سناريوى SSP1-2.6‏ در هر سه دوره آينده 
نسبت به دوره يايه در منطقه موردمطالعه افزايش بيدا خواهد كرد كه ميزان آن به‌طور متوسط در دوره ۲۰۲۱ 


تا ۲۱۰۰ برابر با 71 درصد خواهد بود (شكل ۵). 
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شكل ۵- توزيع فضايى بارش در منطقه شمال غرب كشور بر اساس سناریوی SSP1-2.6‏ در سه دوره آينده 


نسبت به دوره يايه 


توزیع فضایی مقادیر بارش بر اساس سناریوی SSP3-7.0‏ نشان داد كه بارش در منطقه موردمطالعه در سه 
دوره آینده موردبررسی (۵۲۰۲۱-۲۰۵۰ ۲۰۵۱-۲۰۸۰ و ۲۰۸۱-۲۱۰۰) به ترتیب» 0۲۹/۸ ENA‏ و ۸۲/۳ میلی‌متر 
کاهش خواهد يافت كه نسبت به دوره يايه به ترتيب الى 1/۲ و YY/Y‏ درصد کاهش را نشان می‌دهد. 
کمترین تغییرات کاهشی حاصل مربوط به نواحی شمالی منطقه و بیشترین تغییرات مربوط به نواحی جنوب 
منطقه موردمطالعه کاهش بيدا خواهد کرد که اين میزان کاهش در سطح منطقه در دوره ۲۰۲۱ تا ۲۱۰۰ نسبت 


به دوره al‏ به‌طور متوسط برابر با ۵ درصد خواهد بود (شكل A‏ 
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شکل 1- توزیع فضایی بارش در منطقه شمال غرب كشور بر اساس سناريوى SSP3-7.0‏ در سه دوره آينده 


نسبت به دوره پایه 


دوره (۲۰۳۱-۲۰۵۰) بين ۰ نا ۰ ميلى متر خواهد بود كه نسبت به دوره ak‏ ۸ درصد (۱/۱میلی‌متر) 
افزایش خواهد داشت؛ اما در دو دوره دیگر نيز (۲۱۵۱-۲۰۸۰ و ۰--۲۰۸۱) میزان بارش به‌طور متوسط در 
منطقه موردمطالعه به ترتيب» ۲ و ۶/۵ درصد نسبت به دوره ab‏ کاهش خواهد یافت. بیشترین تغییرات 
افزایشی و کاهشی در پهنه موردمطالعه نیز همچون دو سناریوی AS‏ به ترتیب مربوط به نواحی شمالی و 
c £g‏ غربى منطقه خواهد بود. درمجموع ميزان بارش در يهنه موردمطالعه بر اساس اين سناریو در دوره 
۲۱ تا ۰ نسبت به دوره يايه Y/A‏ درصد كاهش پیدا خواهد کرد (شکل (V.‏ 


\M‏ جغرافيا و مخاطرات محيطى شمارة اول 


45°0'0"E 47*0'0"E 49*0'0"E 45°0'0"E 47*0'0"E 49*0'0"E 


1985-2014 2051-2080 


40°0'0"N 
40°0'0"N 


38°0'0"N 


38°0'0"N 


36°0'0"N 
36°0'0"N 


-280 
-400 
-510 
-620 
o 50 100 200 _740 0 50 100 


40^0'0"N 
40°0'0"N 


2021-2050 2081-2100 


38°0'0"N 
38°0'0"N 


36°0'0"N 
36°0'0"N 


0 50 100 200 2 200 300 km 


45?0'0"E 47*0'0"E 49*0'0"E 45*0'0"E 47*0'0"E 49*0'0"E 


شكل ۷- توزيع فضايى بارش در منطقه شمال غرب كشور بر اساس سناريوى SSP5-8.5‏ در سه دوره آينده 


نسبت به دوره يايه 


درمجموع بر اساس نتایج حاصل از بررسى سناريوهاى مختلف بهجز سناریوی SSP1-2.6‏ روند بارش در 
دوره‌های موردبررسى در منطقه موردمطالعه كاهشى است و بيشتر تغييرات به مناطق پربارش منطقه 
موردمطالعه اختصاص دارد. تغییرات افزايشى در برخى مناطق موردمطالعه نيز می‌تواند ناشى از نزديكى به 
منابع آبى و مسير رطوبت درياى سياه و مديترانه و یا كوهستانى بودن و پوشش كياهى و جنكلى آن منطقه 
باشد. مقایسه سناریوهای مختلف نیز نشان می‌دهد که سناریوی SSP3-7.0‏ تغییرات کاهشی بیشتری نسبت 
به So‏ سناریوها OLE‏ می‌دهد و برعکس سناریوی SSPI-2.6‏ تغییرات کمتری را OUS‏ می‌دهد که ناشی از 


ویژگی‌هایی است که هر یک از سناریوها به نمایش می‌گذارند (شکل AA‏ 


سال سيزدهم 


card 


نكرى بارش با استفاده از مدلهاى CMIPG‏ تا بایان قرن ۲۱ در dies‏ غرب ايران 
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۱۹۰ جغرافيا و مخاطرات محيطى شمارة اول 


*- نتیجه گیری 


بيش نكرى تغییرات بارش برای اطلاع از میزان تغییرات آن در آینده و در نظر گرفتن تمهیدات لازم برای 
تعدیل اثرات سوء ناشی از تغیبرات آب‌وهوایی بر منابع آبی و کشاورزی از اهمیت زیادی برخوردار است. 
اهمیت پیش‌نگری بارش زمانی آشکارتر است که در جایگاه‌های جغرافیایی ریزش ناچیز يا قابلتوجه و يا 
ناگهانی حجم زیادی از بارش رخ دهد. لذا با توجه به اهمیت موضوع. در اين پژوهش به ارزیابی ۷ مدل 
GCMs‏ از CMIP6‏ در شبیه‌سازی بارش و همچنین پیش‌نگری تغییرات بارش در منطقه شمال غرب کشور 
تحت سه سناریوی واداشت اجتماعی- اقتصادی (SSPs)‏ در سه دوره زمانی آینده (۰۲۰۲۱-۲۰۵۰ ۲۰۸۰- 
۱ و ۲۰۸۱-۲۱۰۰) نسبت به دوره يايه (۱۹۸۵-۲۰۱۶) پرداخته شد. نتايج حاصل از ارزیابی عملکرد 
مدل‌های GCM‏ موردبررسی با استفاده از شاخص‌های خطاسنجی مختلف نشان داد که مدل MPI-ESM1-2-‏ 
LR‏ در منطقه موردمطالعه از دقت بیشتری نسبت به ساير مدل‌ها جهت شبیه‌سازی بارش برخوردار است؛ 
بنابراین بر با استفاده از داده‌های مدل MPI-ESMI-2-LR‏ و سناریوهای آن اقدام به پیش‌نگری بارش گردید. 
نتایج حاصل از پیش‌نگری تغییرات بارش نيز نشان داد كه بارش بر اساس سناریوی خوش‌بینانه SSP1-2.6‏ تا 
پایان قرن ۲۱ در سطح منطقه حدود ۲/١‏ درصد نسبت به دوره ab‏ افزایش بيدا خواهد کرد. با اين وجود 
به‌جز ols‏ شمالی منطقه موردمطالعه در ساير پهنه‌های موردمطالعه میزان بارش در دوره آینده نسبت به 
دوره پایه کاهش خواهد یافت که بیشترین میزان کاهش آن نيز مربوط به نواحی پر بارش در جنوب غرب 
منطقه موردمطالعه است. به‌جز ایستگاه‌های مریوان» ارومیه و تبریز و زرینه در ساير ایستگاه‌ها افزایش بارش 
معنی‌دار شده است که بیشترین افزايش در سطح معنی‌داری ۰/۱۱ درصد مربوط به ایستگاه اردبیل می‌باشد؛ 
اما بر اساس سناریوهای SSP3-7.0‏ و SSP5-8.5‏ میزان بارش تا OLL‏ قرن ۲۱ به ترتیب ۱۶/۵ و V/V‏ درصد 
کاهش بيدا خواهد کرد و بیشترین کاهش أن نیز مربوط به نواحی جنوب غربی منطقه موردمطالعه است. بر 
اساس سناریوی SSP3-7.0‏ در بیشتر ایستگاه‌های موردمطالعه در منطقه بارش کاهشی و معنی‌دار شده است. 
بیشترین میزان کاهش بارش مربوط به ایستگاه‌های سردشت. خوی و جلفا می‌باشد که در سطح ۰/۰۱ درصد 
معنی‌دار شده است. نتایج حاصل از سناریوی SSPS-8.5‏ تا پایان قرن ۲۱ نيز نشان داد که بارش در منطقه 
به‌جز ایستگاه ميانه در تمامی ایستگاه‌ها دارای روند کاهشی است و میزان اين کاهش نسبت به سناریوی قبل 
اندکی بيشتر است. میزان اين کاهش در مناطق غربی و جنوب غربی منطقه موردمطالعه در ایستگاه‌های 
سردشت. پیرانشهر مهاباده يارس آبادء ماکی جلفا و خوی در سطح ۰/۰۱ درصد معنی‌دار شده است. quU‏ 
حاصل با مطالعات مجدی و همکاران (۲۰۲۲) مبنی بر کاهش بارش در بیشتر مناطق خاورمیانه و شمال 
آفریقا بر اساس متوسط ۲۳ مدل GCM‏ مطالعه حیدری و همکاران(۱۳۹۸) بر مبنای بررسی مدل‌های 


كردش كلى جو و استفاده از سناريوهاى RCP‏ مبنى بر کاهش بارش در دورههاى ) ۲۰۵۱-۷۰ و ۰ 


سال سيزدهم پیش نكرى بارش با استفاده از مدلهاى 081186 تا يايان قرن 7١‏ در dius‏ غرب ايران ۱۹۱ 


۲۰۵۰ و همجنين مطالعات صلاحى و همكاران (۱۳۹۵) مبنى بر کاهش ۲ درصدى بارش در دهه‎ (Ae 
در حوضه آبخیز درياجه اروميه مطابقت دارد. از نتايج حاصل از اين يزوهش مىتوان در مقابله با معضلاتى‎ 
همچون خشکسالی. سیلاب‌های ناگهانی و تخریب‌های زیست‌محیطی و همچنین مدیریت منابع آب؛‎ 
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